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We study the future evolution of quintessence/phantom dominated epoch in modified F{R)- 
gravity which unifies the early-time infiation with late-time acceleration and which is consistent 
with observational tests. Using the reconstruction technique it is demonstrated that there are 
QQ ■ models where any known (Big Rip, II, III or IV Type ) singularity may classically occur. From 

f^ ' another side, in Einstein frame (scalar-tensor description) only IV Type singularity occurs. Near 

C") I the singularity the classical description breaks up, it is demonstrated that quantum effects act 

rvQ , against the singularity and may prevent its appearance. The realistic _F(_R)-gravity which is future 

singularity free is proposed. We point out that additional modification of any _F(7?)-gravity by the 
terms relevant at the early universe is possible, in such a way that future singularity does not occur 
even classically. 
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I. INTRODUCTION 

Q-i| The current interest to modified F(i?)-gravity is caused by its success as the gravitational alternative for dark energy 

^ • (for a review, see p|). The number of studies of dark energy aspects of i^(i?)-gravity [3j S |j| clarified its structure 
and possible cosmological applications. It has been demonstrated [3| that there appears the possibihty to unify the 
description of the inflation and late-time cosmic acceleration within the same theory which has standard newtonian 

\l ■ regime in Solar System. Hence, the whole universe expansion history may be obtained as the cosmological solution 
from some modified i^(i?)-gravity. Recently, the class of viable (non-linear, analytical) i^(-R)-gravities was suggested 
in ref.Q- Such theories and their extensions [a, 01 do pass the local tests and successfully describe the (almost) 
ACDM epoch. Generalizations of these theories proposed in refs.p, M, ll3] keep all their nice viability properties 

^^ ■ but may simultaneously describe the inflation, so that whole universe expansion history: inflation, radiation/matter 
dominance, dark energy epoch follows from the same viable -F'(i?)-gravity consistent with Solar System tests. 

The observational data indicate that current dark universe has the effective equation of state paremeter w being 

^Q ' very close to —1. When w = —1 the universe passes through ACDM epoch. If w is slightly less than —1 then we live 

f^ , in phantom-dominated universe and if w is slightly more than —1 the quintessence dark epoch occurs. In all these 
cases part of (or all) energy conditions are violated. It is known that dark energy universe with the effective phantom 
phase ends up in the future, finite-time singularity |ll| . From another point, the effective quintessence universe may 
end up in more general (soft) singularity. This is also true for modified i^(_R)-gravity modeflj which leads to dark 
?H ] universe with the corresponding (phantom/quintessence) effective w. One can always make the effective phantom 
phase being transient, just changing the classical structure of i^(i?)-gravity by the terms which are relevant at the 
very early universe (see the corresponding investigation in [l2| ). Recently, the appearance of future singularity in the 
specific example of _F(i?)-gravity was re-obtained in refs.[13l.ll4|. 

In the present work we investigate the future evolution of modified gravity which unifies the inflation with late- 
time acceleration. Using the classification of finite- time future singularities proposed in ref. [1^ we describe how all 
types of finite-time singularities may occur in modified gravity. Comparing Jordan and Einstein frames we show 
that the transition of the singularity type occurs: whatever type future singularity appears in Jordan frame ( the 
original i^(i?)-gravity), in the Einstein frame (scalar-tensor description) it shows up as Type IV singularity. The 
classical mechanics analogy for future singularity is given. As classical description breaks up near to singularity, the 
quantum effects should be taken into account while approaching future singularity. It is demonstrated how quantum 
effects prevent the future singularity in modified gravity. The realistic model of non-linear _F'(i?)-gravity which unifies 
the inflation with late-time cosmic acceleration without future singularity and which is consistent with local tests is 
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proposed. The early-time modification of tfie gravitational theory may be always proposed in such a way that future 
singularity does not occur while local tests are not violated. 

II. GENERAL STRUCTURE OF F(i?)-GRAVITY AND TRACE EQUATION 

Let us start from the general (Jordan frame) action of _F(i?)- gravity (for a review, see (l|) with matter: 

^i^(«) = /d'.-V^{^ + /:™} , (1) 

The standard field equations are given by 

^g^.F{R) - R^,F'{R) ~ g^,nF'{R) + V^V,i^'(i?) = -yT(„,)^, . (2) 

Here F(R) is a proper function of the scalar curvature R and Cm is the matter Lagrangian. In ([2|). T(„j)^i, is matter 
energy-momentum tensor. By separating F{R) into the Einstein-Hilbert part and modified part as 

F{R)^R + f{R), (3) 

the trace part of the equation of motion ^ has the following Klein-Gordon equation-like form: 

3n/'(i?) = 1^ ^R + 2f{R) - Rf{R) - k'T . (4) 

The above trace equation can be interpreted as an equation of motion for the non trivial 'scalaron' f'{R) with the 
effective potential V^ff . This means that the curvature itself propagates. The above trace equation is very convenient 
to study the possible instabilities and solar system tests [7;] . Note also that potential which appears in trace equation 
is rather formal, it has no direct physical meaning, like energy density or pressure. 

By introducing the auxiliary field A, one rewrites the action ([1]) of the i<"(i?)-gravity in the following form: 

S = ^ y d^x^^{F'{A) iR-A)+ F{A) + £„,} . (5) 

By the variation over A, one obtains A = R. Substituting A = R. into the action ([5|), one can reproduce the action in 
P]). Furthermore, we rescale the metric in the following way (conformal transformation): 

g,.u~^e''g^, , a = -\nF'{A). (6) 

Hence, the Einstein frame action is obtained: 



Se - ^ld^x^^!^R-^gP'^dpad,a-Via)+Ci^ 



Via) = e^G (e-^) - e^/ (G (e"^)) - ^ - |^ (7) 

Here G (e~'^) is given by solving the equation cr = — In (1 + f'{A)) ~ lnF'{A) a,s A — g (e^'^). 
In the Hu-Sawicki model [a|, when the curvature is large, f{R) behaves as 

f{R) - -2A + ^ . (8) 

Here A, a, and n are positive constants. Then the potential V is given by 

V = 2A + 0{A-'') . (9) 

Therefore, the curvature infinity (R = A —^ oo) surely corresponds to the finite value of the potential. 
We may consider the case that f{R) behaves as 

f{R) ~ /„i?" , (10) 



with positive n and /„ being a constant. The potential is found to be 

n-2 

Then if 2 > n > 0, V becomes infinite when curvature R — A goes to infinity. Therefore, in this case the singularity 
could not be reahzed easily. In the model proposed in [£], f{R) behaves as ([TU)) . Therefore, this case qualitatively 
similar to above behavior. 
In the model Q, 

^(J;^ {R-R,r+'+Rr+' 1 WA (^2. 

/o + /i{(i?-i?„f"+V<+^} /i /o + /i{(i?-i?of"+^+i?r'} ' 

when R is large, f{R) behaves as 

/(^)--^ + 72j;rTT' (13) 

which is almost the same with ([8]) and we obtain the expression of V, which is also similar to ([9]), 

V=^ + M-") . (14) 

/i 

In case ([T51) . l/fi correspond to the effective cosmological constant in the inflation epoch and therefore it is very 
large. Then since V becomes positively large, the singularity could not be easily generated. Hence, from the general 
structure of i^(i?)-gravity in the Einstein frame one may deduce its behavior at large curvature. 

We may include the effect of the matter with a constant equation of state (EoS) parameter w. Since the matter 
density behaves as p = poa^^^^^"^' (po > 0) and p has mass dimension 4, by the scale transformation ([5]), p is 
transformed as 

Here a is the scale factor in the Einstein frame, which is given by a = e'^^^a. Note that, the coordinates effectively 
transform as x'^ — > x^^' = e*^" under the scale transformation, general mass dimension n quantity Q transforms as 
Q ^ Q' ~ er'^'^l'^Q. Since the energy density depends on the scalar field a, in the Einstein frame, the potential could 
be shifted as 

V^V=V^ Cpoa"'^'+"^e-(i-3^«)./2 ^y ^ p^a~^(^+-') {F'{A))^-^^ . (16) 

For the model ([5]), we find F'{A) — ^ 1 then the matter correction does not give a strong effect and the singularity 
could be easily realized. On the other hand, in the model (fTO]) . since F'{R) ~ n/„i?"~^, the potential V blows up as 
V ^ i^(i-3M))(n-i) if y; < x/3 and n > 1. Then the singularity could be prevented if we include the matter. 

III. FINITE-TIME SINGULARITIES IN F(7?)-GRAVITY 
In this section, we investigate i^(i?)-gravity models which generate several known types of finite-time singularities. 

A. Big-Rip type singularity 

As the first example, we consider the case of the Big Rip singularity, where H behaves as 

^ = T^. (17) 

to — t 

Here ho and to are positive constants and H diverges at i = to- In order to find the i^(i?)-gravity which generates 
the Big Rip type singularity, we use the method of the reconstruction, that is, we construct F(R) model realizing any 
given cosmology using technique of ref. [16J. The general F(_R)-gravity action with general matter is given as: 



S = J d^X^/^{F{R) + /:„,atter} • 



(18) 



The action ([18]) can be rewritten by using proper functions P{(j>) and Q{4>) of a scalar field (j)'- 

5= /"d4a:y=^{P((^)i? + Q(0)+£,„attcr} ■ (19) 

Since the scalar field (f> has no kinetic term, one may regard cf) as an auxiliary scalar field. By the variation over </), we 
obtain 

O = F'(0)i? + g'(0) , (20) 

which could be solved with respect to </> as = 4){R). By substituting (/> = 4>{R) into the action ((T9l) . we obtain the 
action of i^(i?)-gravity where 

F{R)^P[<i>{R))R + Q(c^{R)) . (21) 

By the variation of the action (J19p with respect to g^^, the equation of motion follows: 

= - ^9^. {P{c^)R + Q(0)} - R^uP{<i>) + V^ V,P(0) - g^. V2p(0) + ir^, (22) 

In the FRW universe, Eg. ip^ has the following form: 

= -6i72p(0)-Q(0)-6i7^^i|^+p 

at 

= (4iJ + 6i/^)F(^)+g(<^) + 2^^!^+4iJ^^^+p (23) 

By combining the two equations in ((23)l and deleting Q{(j>), we obtain 

dPPichU)) dP(6(t)) ■ , , 

Since one can redefine 4> properly as (j) — (/'('ys), we may choose (/) to be a time coordinate: (j) = t. Then assuming p, p 
could be given by the corresponding sum of matter with a constant EoS parameters Wi and writing the scale factor 
a{t) as a = aoe^*^*' (oq : constant), we obtain the second rank differential equation: 



= 2^^ - 25'(0)^ + 45"(0)F(^) + E (1 + -0 Aoao-^<^+"''^e-^(i+-)^(« . (25) 



If one can solve Eg. ([^5]) . with respect to P{(j)), one can also find the form of Q{4>) by using (P5| as 

g(<^) = -6 ig'm' P(0) - 6ff'(</>)^ + E p.oao ^(^+-)e-3(i+-)'^W) . (26) 

Thus, it follows that any given cosmology can be realized by some specific _F(i?,)-gravity. 

In case of (flT)) . if we neglect the contribution from the matter, the general solution of ((25)l is given by 



F(</>) = F+ (to - <^)- + P_ (to - <t>r- , a± ^ '"" + ^ ^ ^^"° ' ^°' i^ , (27) 

when ho > 5 + 2\/6 or /lo < 5 — 2\/6 and 

P(0) = (to - 0)-('-+^)/^ (icos ((to - ^)i, -^^ + 10^0 -1 ^1 ^^^.^ (^(^^ _ ^)i, -^§ + 10^0-1 ^^ , (28) 

when 5 + 2^6 > ft-o > 5 - 2\/6. Using ((201), dH]), and ^, we find the form of F{R) when i? is large as 

F{R) ex pi^"-/2 , (29) 

for /iQ > 5 + 2V6 or ft-o < 5 — 2-^/6 case and 

F{R) ex pC^o+i)/^ X (oscillating parts) , (30) 



for 5 + 2\/6 > ho > 5~ 2%/6 case. Then V behaves as 

for ho > 5 + 2V6 or ho < 5 — 2^/6 case and 

V -- i?(3-'!o)/4 ^ (oscillating parts) , 



(31) 



(32) 



for 5 + 2\/6 > ho > 5 — 2\f% case. Hence, even if the curvature tends to infinity, the potential V does not always tend 
to infinity. Note also that the potential is often unbounded below. 



B. More general singularities 

Let us investigate more general singularity 

H^hoito-ty^ ■ 



(33) 



Here ho and /? are constants, ho is assumed to be positive and t < to as it should be for the expanding universe. 
Even for non-integer /3 < 0, some derivative of H and therefore the curvature becomes singular. Since the case /3 = 1 
corresponds to the Big Rip, which has been investigated in the last subsection, we assume /3 ^ 1. Furthermore since 
(3 = corresponds to deSitter space, which has no singularity, it is assumed /3 7^ 0. When /9 > 1, the scalar curvature 
R behaves as 



R - 12ij2 - I2hl {to - ty^'^ 
On the other hand, when /3 < 1, the scalar curvature R behaves as 

R^6H <^6ho(3ito~ty'^~'^ 
If we write 



PW 



- p-hoito-<P)-''+'/2{l-f3) 



SW 



Eq. (|25p without matter contribution has the following Schrodinger equation like form: 



^ d^Si<t>) , ^5/3/io,, 
°" d02 +1^ 2 ^'' 


^y^ 


-'-^{to-^ 


When (/) = t — > io, in case /3 > 1, one finds 




2 (^0 0) 


< 


f (io - ^y'' 


On the other hand, in case /3 < 1, we find 




2 (^0 (f>) 


> 


^ito~<t>y'' 



-2/5 



s . 



In any case, Eq. |37)) reduces to the following form: 







Vq (to -^y"s 



when = t — > to. Here 



With further redefinition 



Vo = 

T/ "0 

Vq — -T 



-^ , a = 13+1 when (3 < 1 



a = 2(3 when (3 > 1 
y^(to-0)^-"/^ , 5 = y("/4)(^-"/2)-V, 



(34) 
(35) 
(36) 

(37) 

(38) 

(39) 

(40) 

(41) 
(42) 



"'-a- ^-i il + TlP^'"- («) 



Eq. pOjl has the following form: 

dV I fa^ a\ 1 4Vo 

d^ " I ^Te ~ 4 j ^ + (2 -a)' 

Note that y ^ when (/) -> io if 1 - a/2 > but y -> cx) when -> to if 1 - a/2 < 0. Then if 1 - a/2 > 0, Eq. 
reduces to the following form when — > to^ 

(Pip fa'^ a\ 1 , , 

whose general solution is given by 

ip = Ay"/^-^ + By-°'/^ . (45) 

Here A and B are constants of the integration. On the other hand, if 1 — a/2 < 0, Eq. (|43|) reduces to the following 
form when — > tg: 



<Pip ^ 4K) 
dy^ (a - 2) 



o = 3;| + -^^- (46) 



When Vb > 0, the general solution of (l46l) is given by 



P 



2 /!/ 

y — A cos (wj/) + S sin (wy) , lo = . (47) 

a — 2 

Here A and B are constants of the integration. On the other hand, if Vq < 0, the general solution has the following 
fornr 

y = Ae^y + Be-^y , uj = hE^ . (48) 

a — 2 

From the above analysis, one may get the asymptotic solution for P when (j) ^ t^. 

• /3 > 1 case: From (|¥l]) . a = 2/3 > 2 and therefore 1 — a/2 = 1 — /3 < 0, which corresponds to (|1^ . Since we also 
find Vo > 0, the solution is given by dH]). Then by combining ^, ^, ^, and (gTI), we find the following 
asymptotic expression oi P{(j)): 

(0) ~ e(''«/^(''-i»(*o-0)-^+^(i^_^)/3/2(^ieos(c.(io-0)-^+')+Ssin(a;(to-<^)-^+')) , 

When -> to, ^(0) tends to vanish very rapidly. By using ([201), dH]), and ([261), F(i?) looks like (at large R) 

F{R) oc e^'"'/^^'^"^^HT^J i?-i/4 X (oscillating part) , (50) 

which gives 

^_^-(/.a/2(/3-i))(.r5j^j i?9/4 X (oscillating part) . (51) 

Near the curvature singularity i? — > oo, the potential becomes small exponentially. 

1 > /3 > case: From (|4T|) . we find a = /3 + 1 and therefore 1 — a/2 = 1/2 — /3/2 > 0, which corresponds to 
dm). Since 



^-1 

4 



(-9=/3>0, (52) 



the second term in d45|) dominates when </) ^ to if -B 7^ 0. Then by combining d33|) . d36|) . (|42|) . and d45|l . we find 
the following asymptotic expression of P{(\))'. 

P{^) ^ Be-("»/2(i-«)(*o-0)-^^ (i^ _ 0)('3+i)/8 . (53) 



Therefore, 

F{R) - e-iho/2{l-m(-6l3hoR)^^-^^^'''+^^jf/8 ^ (-54) 



Eg. ((35)) shows that when (j) = t^ to, R-^oo in case /3 > — 1 but R ^ in case /3 < — 1. Therefore we find the 
asymptotic behavior of the potential when R ^ oo as 

V ^ e(''"/2(l-/3))(-6/3/iofl)<'^-i>/(''+i'^9/8 _ ^9/8 ^ (-55^ 

which diverges in the limit of i? — )■ oo. 
• /3 < cose; As in 1 > /? > case, one gets a = (3 + 1 and therefore 1 — a/2 > but since 

--l-(--)=^<0, (56) 

the first term in (|45p dominates when <j) ^ to, ii A ^ 0. Hence, the asymptotic expression of P{(j)) follows: 

P(0) ~ Ae-(''°/2(i-'5))(*o-0)^-'^ (ig _ 0)-('5^-6/3+i)/8 . (57) 

Then F{R) is given by 

FiR) ^ (-6/10/3^?)^^'+'"+')/'^''+'^ g-(ho/2(i-^))(-6/.o,3fl)<^-'/<''+^) _ (58) 

Note that —6hof3R > when /ip: R > 0. If < /3 < —1, R ^ oo when (f) = t ^ to and therefore we find 

^ _ ^2-(/5+l)V8-l/(/3+l)g(fto/2(l-/3))(-6/io/3fi)<''-^>/<''+^> ^ (59) 

which diverges when R ^ oo. On the other hand, if/3< —1, R ^ when 4> ~ t ^ to and the potential has the 
following form: 

y^j^2-iP+ir/8-l/(P+l) ^ (60) 

which could vanish or diverge depending on the value of /3. 



When /3 > 1 in ([33|) . i? behaves as in (p4|) . and when /3 < 1, the scalar curvature R behaves as in (|35|) . Conversely, 
when i? behaves as 

R r^ 6H r^ Rq iP + 1) {to - ty , (61) 

if 7 > 2, which corresponds to f3 — 7/2 > 1, H behaves as 

H^^ito-tr^' ' (62) 

if 2 > 7 > 1, which corresponds tol>/? = 7 — 1>0, iJis given by 

H-^l^^i^o~t)-^^\ (63) 

and if 7 < 1, which corresponds to /3 = j — 1 <0, one obtains 

H^Ho + '^j^^ito-ty^' . (64) 

Here Ho is an arbitrary constant, which is chosen to vanish in p3p . Then since H = a(t)/a{t), if 7 > 2, we find 

ait) ex exp [(^ - 1) \[^ito - tr^'A , (65) 



when 2 > 7 > 1, a(t) behaves as 



a(t) ex exp ( g^^^"_^^ (io - ty^ ] , (66) 



and if 7 < 1 , 



a{t) (X exp ( Hot + ^^ ^^°_ ^^ (to - t)-< ) . (67) 



In any case, there appears a sudden future singularity [17| at i = ip ■ 

Since the second term in ([M)) is smaller than the first term, we may solve ([^5)) asymptotically as follows: 

P-Po{l + ^ito- 0)^-^) , (68) 

with a constant Pqi which gives 

FiR) - Foi? + Fii?2/3/(/3+i) . (69) 

Here i^o and Fi are constant. When > /3 > —1, we find 2/3/ (/3 + 1) < 0, which corresponds to ([8]) by identifying 
n = —2/3/ (/3 + 1) or /3 = —n/{n + 2), which could be confirmed later in ([98]) . On the other hand when /3 < —1, we 
find 2/3/ (/3 + 1) > 2. As we saw in ([29|) . the -F(i?) generates the Big Rip singularity when i? is large. Then even if R 
is small, the F{R) generates a singularity where higher derivatives H diverge. 

Let us also investigate how effective equation of state (EoS) parameter WeS for ((64)) behaves when t ^ to. In the 
Einstein gravity, the FRW equations are given by 

^H' = P, ^^{2H + 3H')^p. (70) 

Since the EoS parameter is defined hy w — p/ p, even for the -F(i?)-gravity, we may introduce the EoS parameter WcS 

by 



Then if /3 = 7/2 > 1, it follows 



9 TJ 

wcff = -l-^^. (71) 



w,s^-l-^{to~tr^+^^-l, (72) 

6tlo 



when t ~* to- If 1 > /? = 7 - 1 > 0, we find 

w,sr^-^ito-t)-'+^^-^. (73) 

oho 

Finally if /3 = 7 — 1 < 0, one gets the expression ((73|) when Hq = 0. When Hq ^ 0, on the other hand, we obtain 

- [to - t] 
'0 



^cff~-l-^(io-t) ''' , (74) 



when Ho vanishes. Then if —1 < /3 < 0, WeS —^ +00 when t — *■ io- On the other hand, if /3 < —1, Woff — * —1. 

Eg. ([TO]) also shows that, even for F(i?)-gravity, we may define the effective energy density pcff and the effective 
pressure PcS by 

p,s = ^H^, p^g = -L(2H + SH^) . (75) 

We now assume H behaves as (p3| . Then if /3 > 1, when t ~> to, a ^ exp(/io (to ^ t) ~ / (/? — 1)) — > 00 and 
Pcff, boffl ^00. li (3 — 1, we find a ^ {to — t) " ^ 00 and poff, |Pcff | ^00. If < /3 < 1, a goes to a constant but 
p, \p\ -^ 00. If — 1 < /3 < 0, we find a and p vanishes but |peff| -^ 00. When /3 < 0, instead of (|551). as in (p5| . one 
may assume 

H^Ho + ho{to~t)-^ . (76) 

Hence, pcff has a finite value 3Hq/k^ in the limit t ^ to when — 1</3<0. If/?<— 1 but /5 is not any integer, a is 
finite and poff and Poff vanishes if Ho = or pefi and poff are finite if Hq ^ but higher derivatives of H diverge. 
In [1^ , there was suggested the classification of the finite-time singularities in the following way: 



• Type I ("Big Rip") : For t ^ tg, a ^ oo, p — > cx) and \p\ — > oo. This also includes the case of p, p being finite 
at L g . 

• Type II ("sudden") : For t ^ tg, a -^ ag, p ^ ps and \p\ — > cx) 

• Type III : For t —> tg, a —f ag, p —i- oo and \p\ — > cxd 

• Type IV : For t ^ tg, a ^ Ug, p ^ 0, \p\ ^ and higher derivatives of H diverge. This also includes the case 
when p (p) or both of them tend to some finite values while higher derivatives of H diverge. 

Here tg, Ug and Ps are constants with ag ^ 0. We now identify tg with Iq. The Type I corresponds to /3 > 1 or /3 = 1 
case, Type II to — 1 < /3 < case, Type III to < /? < 1 case, and Type IV to /3 < —1 but (3 is not any integer 
case. Thus, we have constructed _F(i?)-gravity examples which show any type of above finite-time singularity. This 
is natural because it is known that modified gravity may lead to the effective phantom/quintessence phase 11!] while 
the phantom/quintessence dominated universe may end up with finite-time singularity. 

The reconstruction method also tells that there appear Type I singularity for F{R) — R + aR"' with n > 2 and 
Type III singularity for F{R) = R- /?i?-" with n > 0. 

C. Classical mechanics analogy 

Let us now consider the analogy with the classical mechanics in the description of finite-time singularities in modified 
gravity. 

We start with the trace equation: 

2F{R)^RF'{R)~3nF'{R)=-—T. (77) 

By writing F{R) as F{R) =R + f{R), Eq.^ could be rewritten as 

R + 2f{R)-Rf'{R)~3af'{R)^-^T . (78) 

When the curvature is large, we assume f{R) behaves as ([5]). Then Eg. ([75]) reduces to 

i? + 3an(i?"""i) -0 . (79) 

If 

X ^ R-''-' , (80) 

and the FRW metric with flat spatial part is chosen. Eg. ([75]) has the following form: 

X + 3i/x=^X-'/^"+'^ (81) 

3a 

Note X = corresponds to the curvature singularity _R — > oo. Neglecting the second term in the l.h.s. of ([5T|) . the 
equation looks very similar to the equation of the motion in the classical system with repulsive force (l/3a)x^^/^"+^^ 
First we consider the classical equation of motion: 

.T = ±x-i/("+i) . (82) 

3a 

For Eq. ((82|) . one gets an exact solution: 

x = C (io - t)'^"+'^/^"+'^ . (83) 

Here C and to are constants. Note 2 > 2(7i + l)/(n + 2) > 1. Then x vanishes in a finite time t = to, which corresponds 
to the curvature singularity in (j82p . 

We now investigate the asymptotic solution when the curvature is large, that is, x is small. As there is a curvature 
singularity, one may assume ([55]) . Since 

«._1M___5^, ,84) 

R diverges when (3 > —1 but (3 = 0. 

We now consider three cases: 1) /3 = 1, 2) /3 > 1, 3) < /3 < 1, and 4) > /5 > -1. 
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• In case 1) /3 = 1, since 

and therefore, from ([50]) . we find 
and the l.h.s. of (IST]) behaves as 
but the r.h.s. behaves as 



X ^ (to - tf'"''^ , (86) 

X + 3Hx- {to - tf^ , (87) 

±;^-V(«+l)^(i„_i)-2 ^ (88) 

which is inconsistent since the powers of the both sides do not coincides with each other. Therefore, f3 ^ 1. 

• In case 2) /3 > 1, we find 

R = 12ij2 + 6iJ - 12H^ - (to - ty'^'^ , (89) 

and therefore 

X - (io - t)'''^"+'^ . (90) 

In the l.h.s. of (f8T|) . the second term dominates and the l.h.s. behaves as 

X + mx- iHx - {to - i)'^(^"+i'-i . (91) 

On the other hand, the r.h.s. behaves as 

1 -!/(«+!) ^(i„_t)-2/^ . (92) 

da 

Then by comparing the powers of the both sides, one gets 

/3(2n + 1) - 1 = -2/3 , (93) 

which gives (3 — l/(2n + 3) but this conflicts with the assumption /3 > 1. 

• In case 3) < /3 < 1 or case 4) > /? > — 1, we find 

R = 12ij2 + 6iJ - iJ - (to - t)-'^-^ , (94) 

and therefore 

X~(to-t)^''+^^("+^) . (95) 

Then in the l.h.s. of ([5T|) . the first term dominates and the l.h.s. behaves as 

^ + 3Hx-x-{to-tf^"+'^+"-' . (96) 

On the other hand, the r.h.s. behaves as 

1 -V(«+i)^(i„_t)-/3-i (97) 

6a 

Then by comparing the powers of the left-hand side and the right-hand side, the consistency tells 

P{n + I) + n - I = -f3 - I or f3 = -n/{n + 2) . (98) 

This conflicts with the case 3) < /? < 1 but consistent with the case 4) > /3 > — 1, which could correspond 
to the case in [lji| or to ((55|) . In fact, by substituting ((98| into ((95)) . we get 

X ^ (to - t)2("+i)/("+2) . (99) 

which corresponds to ([55)1 . Since > /? > — 1, this singularity corresponds to Type III in |15l |. 
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Then we found the curvature singularity really appears in the Hu-Sawicki model in a finite time as pointed out in 

We now briefly comment what could happen near the singularity in the Einstein frame ([7|)- Since a ^ and 
therefore e*^ ^ 1 in the scale transformation ([5]), we may identify the time coordinate in the Einstein frame and the 
time coordinate in the original F{R) frame. Since ([6]) shows also 



^n+l Jln+1 



by using (|M|) and (|98p. one obtains 



and therefore 



(6/1, 



^y^TT (^0 - ^) 



2(n+l)/(n+2) 



(100) 



(101) 



a 



,^ '"("+/)",, (to - t)"/("+^) . (102) 

(6/io)"+H" + 2) 



On the other hand, the potential V{(t) in ([7]) has the following form: 

Via) - 2A - ^-!^ ^ 2A - (^ (to - t)^"/^"-' . (103) 

Then in the Einstein frame, the energy density pa- and pressure Po-, which are given by 

Pa = \o^ + V{0) , Va = ^^' - V{o) , (104) 

and therefore the curvature in the Einstein frame Re = — ('*^/2) T = (k^/2) (p — 3p) are finite even if t ^ to since 
0-2 ^ y _ 2A ~ (to — t) "''■"+ ■' vanishes. As the exponent 2n/in + 2) is fractional in general, however, the higher 
derivative of H could diverge. Then in the Einstein frame, the finite-time singularity appears as Type IV singularity. 
Thus, we discovered the possibility of change of finite-time singularity nature with the change of frame in modified 
gravity. This is not strange. Indeed, it was shown some time ago 18] that Big Rip singularity which appears in scalar- 
tensor theory qualitatively changes its form in the equivalent Jordan frame. Exactly in the time when scalar-tensor 
theory becomes singular, the mathematically-equivalent _F(i?)-gravity (which is not completely physically equivalent 
to scalar-tensor theory ,19]) becomes complex |la |. 

The important remark is in order. The general (Jordan frame) action of i^(i?)-gravity with matter in (HI) can 
be rewritten in the Einstein frame action ([7]). One should note that the sign in front of the kinetic term in ([7]) 
is always canonical, and therefore there never appears phantom in the transformed Einstein frame. Note that this 
transformation is exact one. 

In the Einstein frame, the FRW equations have the following form: 

\H^=Pa + p, ~^(3H^ + 2h) =p,+p. (105) 

Here p and p express the contribution from the matter and pa and pa are defined by (|104p . If we neglect the 
contribution from the matter, the EoS parameter is given hy w = Pa/pa- and therefore w > — 1 if we assume V{a) > 0, 
which corresponds to the fact that a cannot be phantom in the Einstein frame. 
When the Hubble rate H is given by 

H^ho (to - ty" , (106) 

the effective EoS parameter looks like 

Since we are considering the period t < to and the universe is expanding, we assume /lo > 0. Then WcS can be greater 
than —1 when /3 < 0, which corresponds to the Type II (0 > /3 > —1) or Type IV (/3 < —1). Since the EoS parameter 
is greater than —1 in the Einstein frame. Type I (Big Rip type /3 > 1) and III (1 > /3 > 0) singularities are not allowed 
(the phantom absence) but there can occur Type II or IV singularity. 
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Even if there is no singularity in Einstein frame, there could appear a finite time singularity in the Jordan frame. 
We may consider the following action in the Einstein frame ([7]), where the potential is given by 



V{a) = Vbe V ^1-0 . (108) 

Here Vq and Hq are positive constants. Then by assuming the FRW universe in the Einstein frame 

dsl ^ -dtl + aE {tEf Yl (^^')^ ' (109) 

1=1,2,3 

one obtains a solution: 



(aE(tE)oci^") , a=y2^1n^. (110) 



hn 

TT ^ 

<E \ ^ ^ "' '^ J ' V 3 ^EO 



Then there is no singularity although a becomes infinite in the limit of Ie —* +oo. Since the metric tensor gEi_Lv in 
the Einstein frame is related with the tensor gj^^^ in the Jordan frame by Qeiiu = c^'^gjp.i^, the metric in the Jordan 
frame is given by 

dsl=e-'' i-dtl, + aE{tEf Y. (^^')' • (111) 

\ 1=1,2,3 / 



Since the FRW universe in the Jordan frame is given by 

1=1,2,3 

by comparing (jllip with (J112p . the time coordinate tj in the Jordan frame is given by 



ds^j ^ -dt^j + aj (tjf J2 {dx'f ^ (112) 



,-...±/.ee-".-..^(|L)- 



(113) 



6 

One now choose the sign ± so that the directions of the time in the Einstein frame and the Jordan frame could not 

ho. 
6 



be changed. If 1 — -^ < 0, i.e., ho > 6, the limit t — > oo corresponds to tj -^ tjQ. Then we find 



aj(ij) = e--/^aE(iE)ociE""-«(tjo-tj)^-^ . (114) 

Hence, if ho > 6, aj (ij) diverges at tj — > ijo, that is there appears a finite singularity even if there is no singularity 
in the Einstein frame. 

D. Quantum effects near finite-time singularity 

Near the future singularity at t = to) the curvature becomes large in general. As a result, near the singularity 
the quantum fields/quantum gravity effects become very important again. All classical considerations are not valid, 
all speculations about future cosmic doomsday cannot restrict the classical theory structure because quantum ef- 
fects can stop (or shift) the future singularity |20|. Moreover, the quantum corrections usually contain the curvature 
powers, which become important near the singularity. Hence, any claim about the appearance of the effective phan- 
tom/quintessence phase in modified gravity which subsequently enters the future singularity is not justified without 
quantum effects account near to singularity. One may include the massless quantum effects by taking into account 
the conformal anomaly contribution as back-reaction near the singularity |20|. The conformal anomaly Ta has the 
following well-known form: 

Ta = 6 ( F -f |ni? j + 6'G + 6"ni? , (115) 
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where F is the square of 4d Weyl tensor, G is Gauss-Bonnet invariant, which are given by 

r — -it — Zn^vii + -K^iypcr-K , 

G = R^- AR^.R'"' + R^.upaR'"''"' ■ (116) 

In general, with A^ scalar, A^i/2 spinor, iVi vector fields, N2 {— or 1) gravitons and -/Vhd higher derivative conformal 
scalars, h and h' are given by 

AT + 6iVi/2 + 12iVi + &nN2 - 87Vhd 
^" 120(47r)2 ' 

JV + lliVi/2 + 62A^i + 1411iV2 - 28A^HD 

360(47r)2 ■ ^ ' 

As is seen 6 > and 6' < for the usual matter except the higher derivative conformal scalars. Notice that h" can be 
shifted by the finite renormalization of the local counterterm R^ , so h" can be an arbitrary coefficient. 
By including the trace anomaly. Eg. ((78)) is modified as 



R + 2/(i?) - RJ'iR) - 3n/'(i?) = -^ (T + Ta) . (118) 

For FRW universe, we find 

F = , G = 24 (hh'^ + H^\ . (119) 



We now assume H behaves as (|33p and neglect the contribution from matter by putting T = 0. Then in case of the 
Hu-Sawicki model, we assume f{R) behaves as ([5]) and Eq. (|118p becomes 



First we consider the case that 26/3 + 6" = and therefore Ta = G. li —1 < f3 < 0, R behaves as (to — t)~ and G 
behaves as G ~ 2AHH^ ^ {to — t)^ . Since —3/3 — 1 > — /3 — 1 when /? < 0, Ta can be neglected compared with 
R. Here, D (an/i?"+^) ~ (to — ty '^ ' . Then by comparing the r.h.s. with the l.h.s. in (|120p . one reobtains 
(j98p . Therefore, the curvature singularity appears in a finite time and the quantum correction does not prevent the 
singularity when 26/3 + 6" = 0. We may, however, consider the case 26/3 + 6" 7^ 0. In this case, Ta behaves as 

Ta - 06 + h''\ UR ^ (to - t)^^-' . (121) 

Since R ^ {to - ty'^^'^ and D {an/R^+^) - (to - t)(/3+i)(«+i)-2^ ^-^^ terms in the l.h.s of ^20]) are always less singular 
than Ta since 

-P-l, (/? + 1) (71 + 1) - 2 > -/? - 3 . (122) 

This indicates that Eq. (|120p docs not allow the singular solution and the curvature singularity does not appear. 
Therefore, in case 26/3 + 6" 7^ 0, the quantum effects prevent the singularity appearance. 

In the above analysis, DR term acts against the singularity. The OR term is generated by a local term R^, which 
shows that if one modifies F{R) or f{R) by adding R^ term as 

F{R)~^ F{R)+-fR'^ , (123) 

with a constant 7, the curvature singularity could not be generated. Note also that near to singularity, quantum 
gravity effects become dominant. Then, that is quantum gravity which is not constructed yet should define the 
universe behaviour near to singularity. 

Actually, this observation was made some time ago in ref. [12| where it was shown that in the presence of R^ term 
(or higher powers of curvature) as well as negative power curvature terms the effective phantom phase is transient 
and future singularity never occurs. This suggests the scenario to avoid the future singularity by the modification of 
early-time relevant part of i^(-R)-theory, for instance by powers of curvature. Of course, such modification should not 
destroy the inflationary phase and its exit as well as perturbations structure. It is also important that by construction 
such modification does not violate the local tests viability of modified gravity. 
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IV. REALISTIC MODEL WITHOUT FUTURE SINGULARITY 

It has been shown that if f{R) behaves as ([8]) for large curvature there could occur the curvature singularity in a 
finite time. In the Hu-Sawicki model, the behavior of f{R) in ^ generates the late-time acceleration. 

Note, however, the late-time acceleration occurs if /(i?o) becomes almost constant in the present universe: 

/(i?o) = -2i?o , f'{Ro) ~ . (124) 

Here Rq is the current curvature of the universe and we assume Rq > Rq. Almost constant f{R) plays the role of the 
effective small cosmological constant: A; ^ — /(i?o) = 2i?o, which could generate the accelerating expansion in the 
present universe. 

In the Hu-Sawicki model [5| the following condition is also satisfied: 

]imf{R} = 0. (125) 

The condition shows that there occurs a flat spacetime solution, where R = 0. 

Then we need not to require the behavior ([8|) when the curvature is really large R ^ Rq- For example, if f{R) 
behaves as 

/(i?) = -(l-77)i?-2r;Aj + 0(i?-i) , (126) 

we find 

F{R)^r]{R-2Ai) + 0{R-^) . (127) 

In (|126[) . 77 is a positive constant. Eq. (|127p tells that the effective gravitational coupling Koff is given by k^jj — tjk^ 
and the effective cosmological constant, which could generate the inflation, is given by A/. For ()127p . the potential V 
in 17]) is given by 

V^^. (128) 

77 

Then if A/ is large enough, the curvature singularity could not be easily realized. Hence, if F{R) goes to the Einstein 
action (the gravitational constant could be changed) with the large cosmological constant, the value of the potential 
for large R becomes the order of the cosmological constant. Then if the cosmological constant is large enough, the 
value of the potential becomes large and the value cannot be easily realized. 

An example satisfying the conditions (|124p . (|125p . and (|126p could be given by (compare with general viable 
construction [W^ ) 

}{R) = fi{R) + fL{R), 

fj{R) ^ _ (1 _ ,,) i? _ 2ryA/ + 2rjAje^'-'^W^^^r-{i+m-vfRV8v'A-, ^ ^^(^^ ^ -2ai?2e-^"'«' . (129) 

Here c, a, and ^ are positive constants satisfying the conditions A/ ^ 1/a and c ^ 1. When the curvature has 
the order of 1/a, fi{R) behaves as fi{R) ~ R x O {R/Aj) and therefore fi{R) could be neglected compared with 
fL{R) ^ R- Then since 

f'{R) - fLiR) = -4ai? (1 - ca^R^) e'^^'^' , (130) 

compared with (|124p . we find 

i?o - -^ , Ro-^ — e-^ . (131) 

y/ca ca 

Since about 70% of the energy density in the universe is dark energy, we find Rq ~ 0.7i?o- Then one has 

c- (0.7)^6-2-0.07 . (132) 

We should note that the model (|129p unifies the late acceleration and the inflation in the early universe like models 
[8|, Isl and the future singularity does not appear. 
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We now check the model (J129I) could satisfy the constraints on the correction to the Newton law. It is convenient 
to work in the Einstein frame ([7]). The potential V{a) in ([7]) gives a mass ttt-o- for the scalar field a: 

, _ 1 d^V{a) _lj A AFjA) 1 | 

'^"2 da2 2\f'{A) (F'{A)f F"{A)j ^ ^ 

In order that the correction to the Newton law could be small, the mass m„ should be large enough. The Newton 
law has been checked on the earth or in the solar system. In air on the earth, the scalar curvature could be given by 
A = i? -- 10~^° eV^. On the other hand, in the solar system, we find A = i? -■ lO"*^^ eV^. Since Rq ^ lO"*^'' eV^, one 
may assume A/ 3> i? 3> I /a. Then 

and therefore 

l//(i?)l , l/L(i?)l « i? , i/;(i?)i,i/i(i?)i«i. (135) 

Note that f'/{R) is negative and if A/ - (lOi^-^^ GeV)^ = 10'**^-54eV^ we have 

/l'(i?)^ (10^*^-5^ eV')~' . (136) 

We also find 

which is negative. Then in the air of the earth, we find 

/i'(i?)^(lOi°"ev)"\ (138) 

which is extremely small and in the solar system, 

/i'(i?)~(lOi""ev)"' , (139) 

which is also extremely small. Then the mass of scalar field a could be given by 



m 



2., ^ ., ^ ., 1048-54 .r2 _ /inl5-18 



2F"{R) 2f'/{R) 



1048-54gy2^^^Ql5-18GeV)\ (140) 



which is very large and the correction to the Newton law is beoynd of the observational capacities. Hence, the viable 
modified gravity which unifies the inflation with late-time acceleration and which belongs to general class of ref. [lO| 
is free of future singularity. 

As it was shown in the arguments around (fTO]) . pTJ) or (|123p . instead of (|126|) . if we add the term i?" (0 < n < 2), 
the curvature singularity could be avoided in the same way as in |12l |. 

Eq. (fTT|) tells that there could be finite future singularity. We also find that the Big Rip singularity could be realized 
in F{R) ^ i?" {n > 2) theory as shown in (|^^ . This could be confirmed by the trace equation (JH). We now assume 
(TTOl) with n > 2 and therefore F{R) ~ /ni?" with n > 2 when the curvature is large. Then the trace equation ((4| 
reduced in the following form: 

The explicit solution of (|14ip is given by a Big Rip solution: 

h^!!'^ , , 3n2 - 4n + 2 ± J (3n^ - 4n + 2? + 8(n -2)(n- l)n(2n - 1) 

H = ^^^- , 4") = ^ ^, ', ^ — . (142) 

ta-t ^ 4(n - 2) ^ ' 

Since 3^^ — 4n + 2 = 3(n — 2)^ + 8(n — 2) + 4, h]^' is surely real and positive. Then there could appear finite future 
Big Rip singularity in F{R) ^ R" (n > 2) theory. However, considering several such terms with different powers of 
curvature may often lead to stable de Sitter solution which prevents the evolution to singularity. 
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V. DISCUSSION 

In summary, we studied the future evolution in the models of i^(i?)-gravity consistent with local tests and unifying 
the universe expansion history. Working in both (Einstein and Jordan) frames we demonstrated how Big Rip or type 
II, III and IV future singularities appear. Using the reconstruction method [I^ it is shown which models may lead to 
any specific type future singularity. It is remarkable that whatever type singularity may occur in Jordan frame, in the 
corresponding Einstein frame the only IV type singularity appears. It is also interesting that as in the scalar-tensor 
theory, the appearing singularity does not lead to the violation of local tests when modified gravity is consistent 
with them. Nevertheless, some indications to possible future singularity may be expected to be found in the current 
observational boundspj|.l2l|. 

Even if _F(i?)-gravity under investigation develops the future singularity, there is always the way to modify its 
structure by terms relevant at the early universe. Such terms, like in the model [12|, prevent the development of the 
singularity. The explicit example of realistic _F'(i?)-gravity which does not lead to future singularity even classically 
is proposed. ^From another point, the classical description breaks down near to singularity. The account of quantum 
effects becomes necessary in this situation. Taking into account the quantum effects of conformal fields, we show that 
such quantum effects naturally act against the singularity appearence. Of course, in order to understand how realistic 
is our quantum description of future singularity one should take into account quantum gravity. 

As the final remark, let us note that other modified gravities (Gauss-Bonnet gravity, string-inspired gravity, theory 
with inhomogeneous equation of state) [l| also naturally describe the quintessence/phantom dominated universe which 
may evolve to future singularity. It would be of interest to study the structure of future singularities in such models 
too. 
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